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FOREWORD 

This report documents the results of a study on high frequency 
cutoff filter analysis and was performed by Lockheed's Huntsville 
Research & Engineering Center for NASA-Marshall Space Flight 
Center under Contract NAS8-30515. 

The present study was performed during the period 1 July 
1971 to 6 March 1972 in accordance with Modification 9 of Appen- 
dix F. A second independent study — "Payload Deployment Study," 
was completed under this same contract and is documented in 
Volume I. 

The NASA- MS FC technical coordinators for the study were 
Mr. S. N. Carroll and Mr. H. H. Hosenthien of S&E-ASTR-A. 
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SUMMARY 

The problem of digital implementation of a cutoff filter developed by 
Lockheed Missiles and Space Company's Huntsville Research and Engineering 
Center is approached with consideration to word length, sampling rate, accuracy 
requirements, computing time and hardware restrictions. 

Computing time and hardware requirements for four possible pro- 
gramming forms for the linear portions of the filter are determined. Upper 
bounds for the steady state system output error due to quantization for digital 
control systems containing a digital network programmed both in the direct 
form and in the canonical form are derived. This is accomplished by 
defining a set of error equations in the z domain and then applying the final 
value theorem to the solution. Quantization error was found to depend upon 
the digital word length, sampling rate, and system time constants. The error 
bound developed may be used to estimate the digital word length and sampling 
rate required to achieve a given system specification. From the quantization 
error accumulation, computing time and hardware point of view, and the fact 
that complex poles and zeros must be realized, the canonical form of pro- 
gramming seems preferable. 

Finally, the basic input-output characteristics of the digitized high 
frequency cutoff filter is investigated. The filter has demonstrated its ability 
to pass signals with minimum distortion in a given frequency band while mini- 
mizing signals whose frequencies lie outside the band. 

For dc or dc plus a single high frequency signal, the dc signal is 
transmitted and the high frequency signal is blocked by the digital filter. When 
a low frequency signal, or a low frequency plus a dc signal appears at the input, 
the passage is characterized by hysteresis. This is not critical and is shown to 
be insignificant at frequencies less than 0.4 of the cutoff frequency. 
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The phase lag due to hysteresis for a low frequency signal decreases 
with decreasing input frequency, f , below cutoff, f c (from 11.5 for f = 0. 9 f c 
to 0.7° for f = 0.2 f c ), while the output-input amplitude ratio increases with 
decreasing input frequency below cutoff (from 0.26 for f = 0.9 f c to 0.986 for 
f = 0.2 f c ). 

If the input contains multiple high frequency components in addition to 
a low frequency component, the latter is passed in a "qua si -sampled" manner. 
When a high frequency signal alone appears at the input, sharp cutoff occurs 
and the input signal is subdued. If the input consists of two out-of-phase high 
frequency signals, a periodic step waveform appears and the magnitude of the 
output signal is very small compared to the magnitude of the input waveform. 
The period of the output waveform is equal to the least common multiple of 
the periods of the high frequency signals and the sampling period. For effi- 
cient filter operation, studies show that the cutoff frequency should be in the 
neighborhood of one half the lowest expected frequency above cutoff. 
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Section 1 
INTRODUCTION 

In most aerospace control problems, it is desirable to have a filter 
which passes signals with minimum distortion in a given frequency band 
while minimizing signals whose frequencies lie outside the given band. One 
simple filter which has performed this task quite well in a number of simu- 
lations is the high frequency cutoff filter containing a diode bridge developed 
by Lockheed- Huntsville. In previous studies (Refs. 1, 2 and 3), Lockheed had 
simulated the filter in both analog and digital forms. The resulting Simula 
tion methods were primarily intended to be used in evaluating control loop 
performance on Saturn vehicles in the presence of winds. These studies 
have shown the cutoff filter to be as good or better than linear filters which 
were designed around specific frequencies, while the cutoff filter is virtually 
independent of the frequency of disturbance (so long as it lies sufficiently 
above the cutoff frequency). The primary purpose of this report is to develop 
a basic analysis technique for evaluating equipment requirements in order to 
obtain satisfactory performance, and to investigate the basic characteristics 
of the digitized high frequency cutoff filter. 

The main body of this study is divided into four sections. The layout 
of subsequent material is as follows: Section 2 describes an analog version 
of the high frequency cutoff filter (HFCF), develops its digital implementation, 
and presents a technique for evaluating equipment requirements. Section 3 
presents the basic input-output characteristics of the HFCF. Finally, 

Section 4 embodies the conclusions and recommendations. 
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Section 2 

DIGITAL IMPLEMENTATION OF THE HIGH 
FREQUENCY CUTOFF FILTER 

2.1 BASIC DESCRIPTION OF THE HIGH FREQUENCY CUTOFF FILTER 

A block diagram of the continuous high frequency cutoff filter (HFCF) 
is shown in Fig. 2-1. The purpose of the HFCF is to minimize all frequency 
components in the incoming signal which are above a certain value, C0 c> 
Signal components below the cutoff frequency, 0) , are to be passed with a 
minimum of distortion. To accomplish this, the filter utilizes a diode bridge 
in conjunction with the linear networks, as illustrated in the figure. The in- 
coming signal, e in , is first passed through a second-order lag network, then 
split into two paths. The purpose of this splitting is to develop two signals, 
e. and e 0 , which are in phase at low frequencies (below CO ) and 180 degrees 
out of phase at high frequencies (above 0 oj. When the high frequencies are 
passed to the diode bridge, e^ appears at when e^ is positive and at x^, 
when ej is negative. This minimizes passage of the high frequency signals. 
When low frequency signals occur at the bridge, the signals pass with a min- 
imum of distortion. 

2.2 DIGITAL CONFIGURATION 

A block diagram of the digital configuration of the continuous HFCF is 

shown in Fig. 2-2. The incoming continuous signal, e in (t), is converted by 

the analog to digital converter to a sequence of numerical values denoted by 

e. (n), where 
in 


e. (n) = e. (nT), n = 0, 1,2, 
m m 


(1) 
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Filter Leg 1 




CD, < 2 CD 1.5C0 < ox, < 3oo 

1— c c— 2— c 

l 1 < 0.6 3.0<?, 2 <0.8 


e o (s) 


Fig. 2-1 - Block Diagram of Continuous High Frequency Cutoff Filter 
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The operators D^(z) and D^z) symbolically represent linear recursion 
formulas or difference equations that can be programmed in a computer. 

These formulas perform the same signal- splitting and phase- shifting func- 
tions as the linear networks in the analog configuration. The block marked 
digital logic represents a set of logic statements that will perform the same 
function as the diode bridge and output amplifier for the analog filter. 

The method used to obtain the difference equations and the development 
of the logic are discussed below under individual headings. 

2.3 SYNTHESIS OF LINEAR DIGITAL FILTER 

One useful approach in designing a linear digital filter involves finding a 
set of difference equations having a system impulse transfer function that 
closely resembles a known analog filter transfer function. The method used 
for determining the impulse transfer functions represented by D^(z) and D^z) 
in Fig. 2-2 is based on a substitution rule developed by Shaw (Ref. 4). The subs 
titution rule was derived by comparing the classical results of linear differen- 
tial and difference equations. Thus, the transient responses of filters related 
by the substitution rule agree exactly at the sampling points. Consequently, 
the digital filter corresponding to a stable analog filter is stable. 

The transition from differential equation to difference equation (Ref. 4) 
is made in the following manner. 

Substitution Rule 


Consider the differential equation 

f(D) y(t) =Kg(D)x(t) (2) 

where D is the differential operator, Dy(t) = dy/dt, f and g are polynomials 
of degree p and p\ respectively, each with leading coefficient unity. The 
functions x(t) and y(t) are filter input and output, and K is a gain constant. 
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Next, consider the difference equation 


f (E) y(n) = K g(E) x(n) (3) 

where E is the shifting operator, Ey(nT) = y(nT + 1), and T is the sampling 
interval. The functions f and g are polynomials in the operator E. The 
steps in the transition from Eq. (2) to Eq. (3) are as follows: 

(1) Express f{D) and g(D) as products of linear factors, i.e., 

P P' 

f(D) = 77 (D-p.), g(D) = 77 (D-p) 

i=l j=l J 

(2) Replace D by E 

Tp 

(3) Replace each root p by e r 

1/2 p " p ’ 

(4) Replace K by K(TE A/ ^) 

(5) Expand the result, i.e., 


— P TP- 

f<E) = 77 (E- e 


— P' TP- 

g(E) = 77 (E- e 

j 


K =K(TE 1 ^ 2 ) 


P - P' 


The multiplier of K, E^ p p " ‘ 2 , is the proper phase adjustment which adds 
(p-p')/2 to all input subscripts and amounts to central approximation. 
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Impulse Transfer Function 

BasOd upon the substitution rule, the impulse transfer function, D(z), of 
the digital filter is obtained from the transfer function, F(s), of a known analog 
filter via the correspondence 


F(s) = K 


P' 

i?i (s ~ Pi> 

P 

n (s-p.) 

i = 1 


K{Tz A/ ) 


n Tp i 

1/2 P-P' i = l {z ' e > 


P Tp 

77 (z - e x ) 
i = 1 


= D(z) (4) 


Example: 


The transfer function of the analog filter Leg 2 is 


F(s) = 


1 + s 2 /co 2 


(s /oij + s/oo 1 + 1) (s^/co^ + 2 ^ S/co 2 + 1) 


2 / 2 


with 


2 2 2 

(w 1 u 2 /co c ) (s + jo> c ) (s - jco c ) 
(s - pj) (s - p 2 ) (s -p 3 ) (s - P 4 ) 


P\ - ~ + j 

p 2 = " ^“l ' j ( “i 

P 3 ^2 W 2 + j ^ W 2 
P 4 = ^2^2 " j ^°2 


2 _ 2, 

-o^i) 


1/2 


1/2 


2 

’“I 

-> i/ < 

2 

,2 ■ 1/2 

-“2 

^2 

2 , 

, 2> 1/2 

" w 2 

» 2 > 


(5) 
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where co is the cutoff frequency. The correspondence Eq. (4) leads to the 
impulse transfer function 


D{z) = 


(03 2 a} 2 /co 2 ) (T z 1//2 ) (z - e c ) (z 


e ) 


Tp 


T P 3 


(z - e ) (z - e ) (z 


Tp, 


i T P 4 

)(z-e 4 ) 


( 6 ) 


Expanding Eq. (6) yields 


or 


D(z) 


.2 2/2. ,_2. . 3 

(oo 1 (j0 2 /a) c ) (T M^z 

471 3"77 2“ 

z + djZ + d 2 z 


+ c^ z + c^z) 

+ d- z + d . 

3 4 


D(z) 


(oo^M 2 /oo^) (T 2 ) (c 1 z' 1 + c 2 z 2 +c 3 z 3 ) 
1 + djZ -1 + d 2 z“ 2 + d 3 z" 3 + d 4 z“ 4 


( 7 ) 


where 


1 

-2 cos a) T 
c 

1 


“^2^2^ l/2 ”^l^l'^' l/2 

-2 (e cosoo 2 (1 - ; 2 ) ' T+ € 1 1 coscOjU-^) 1 '^) 

-2t, ,0) T -2t, ? (0_T -(£,.0) +t, ? 0) ? )T ,/ 2 

e 11 +e 22 +4e 11 22 (ccsa^l - ^) 1/2 T) 

(cosco 2 (1-; 2 ) 1//2 T) 
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“ .03. + C, 2 C0 2 ) T l/2 

d 3 = -2(e 11 L L cosu 2 (1 -i 2 ) ' T + 


(2; 2 w 2 + ^ 1 co 1 )T 


l/2 

cosaij (1 -4 1 ) 1/ T) 


d. 


- 2 <tA + i 2 “2 )t 


In this study, the digital filter coefficients were normalized to have a 
gain of one at zero frequency. The coefficients of the linear portion of the 
HFCF are tabulated in Table 2-1 for the following parameters. 

03 =3.0 rad/ sec 

c 

(Oj = 4.5 rad/sec, =0.5 

( 0 2 = 6.0 rad/sec, £, 2 =0.6. 

2.4 FREQUENCY RESPONSE OF LINEAR FILTER 

Replacing s by j03 in Eq. (4) gives the frequency- characteristic function 
of the analog filter as F(jco). Hence, the log magnitude and phase functions 
are 20 log | F(jco ) j and tan * [im F(jo3 )/Re F(jco )} , respectively. 

Similarly, replacing z by e in Eq. (4) gives the frequency character- 
istic function of the digital filter as D(e The corresponding log magnitude 

and phase functions are 20 log |D(e^ WT )| and tan 1 [im D (e ^ T )/Re D (e ^ WT )] , 
respectively. 

The frequency characteristics of the linear networks of the high frequency 
cutoff filter are plotted in Figs. 2-3 through 2-6. The log amplitude curves for 
the linear digital networks are very close to that of the analog networks except 
near the Nyquist frequency (i.e., half of the sampling frequency). Such departure 
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Table 2-1 

LINEAR FILTER COEFFICIENTS 


Sampling Interval T = .04 sec 

■BBS 



Filter 
Leg 1 

Filter 
Leg 2 

Filter 
Leg 1 

Filter 
Leg 2 

Filter 
Leg 1 

Filter 
Leg 2 

a 

o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

a l 

0.0 

0.1024040 

0.0 

0.3676234 

0.0 

1.131135 

a 2 

0.0014728 

-0.2033352 

0.0250528 

-0.7101939 

0.6626039 

-1.599667 

a 3 

0.0 

0.1024040 

0.0 

0.3676234 

0.0 

1.131135 


-3.505654 

-2.883835 

-0.8213607 

b 2 

4.654637 

3.276310 

0.5995070 

b 3 

-2.773765 

1.727651 

0.1622874 

b 4 

0.6262534 

0.3602284. 

0.0467448 


Form of Impulse Transfer Function: 


D(z) = 


x -lx -2 X " 3 
a + a . z +a 0 z +a-,z 

o 1 2 3 

- 1 -2 -3 -4 

1 + b^z + b£Z + b^z + b^z 
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Fig. 2-3 - Magnitude Comparison for Filter Leg 1 
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Fig. 2-4 - Magnitude Comparison for Filter Leg 2 


12 


LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 


Phase (deg) 


LMSC-HREC D2 2 5 642 -II 



Fig. 2-5 - Phase Comparison for Filter Leg 1 
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of the log magnitude curves is inherent in the sampling system since the slope 
of the log magnitude function must be zero at the Nyquist frequency. 

2.5 PROGRAMMING METHOD (References 5, 6 and 7) 

Four possible programming forms were considered for implementing 
a pulse transfer function D(z): 

1. Direct Programming 

2. Canonical Programming 

3. Cascade Programming 

4. Normal (or Parallel) Programming. 

To determine a choice between these forms, some basis of comparison must 
be used. The scheme used for comparison in this study is based on the numbe 
of registers (memory size) required for data and coefficient storage, the com- 
puting time involved in additions, multiplications and data transfers, and error 
accumulated due to input, output and arithmetic quantization. 

2.5.1 Direct Programming 

An impulse transfer function D(z) containing n poles and m zeros can 
be expressed in the following general form 


D(z) = 


E ^l z ) a n + a.z ^a^ 2 +... + a z 1X1 
o _ U 1 2 m 

E i 1 + b.z" 1 + b 7 z' 2 + b z" n 

12 n 


, m < n 


(8) 


The difference equation which describes the direct programming form can be 
obtained by cross multiplying Eq. (8) and taking the inverse transform. Thus 


m 


: o (k) = £ a j e i( k -j)-Z b i e o< k -j> (9) 

j =0 j = 1 
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where e^(k) denotes the input evaluated at time kT, and T is the sampling 
interval. A general block diagram of this form is shown in Fig. 2-7. If all 
the a^ and b, are not zero and one, the number of storage registers for direct 
programming is 


Data 

m 

+ 

n 

Coefficients 

m 

+ 

n+ 1 

Multiplications 

m 

+ 

n+ 1 

Additions 

m 

+ 

n 

T ransfers 

m 

+ 

n 


2.5.2 Canonical Programming 

To derive the difference equations which describe this form, one may 
introduce an intermediate variable F(z) into Eq. (8). Let 


E o (z) £ 


F(z) 


j=0 


a. z 
J 


-J 


F(z) 

E^z) 



The inverse z-transform of Eqs. (10) and (11) are then 

m 

e (k) = S a ; f ( k " j) 

j =0 J 


f(k) = e.(k) - ^2 b.f( k ~j) 

j = l 3 


( 10 ) 


( 11 ) 


( 12 ) 


( 13 ) 
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A general block diagram of this form 

the a. and b. terms are not zero and 
J J 

for canonical programming is 

Data 

Coefficients 
Multiplications 
Additions 
T ransfers 


is shown in Fig. 2-8 for m =n. If all 
3ne, the number of storage registers 

n 

m + n+ 1 
m + n + 1 
m + n 
n 


2.5.3 Cascade Programming 

The difference which describes the cascade programming form can be 
obtained by factoring Eq. (8) into the following form for m = n : 


D(z) = 


!+ e 2 z 


-1 


1 + e z 
n 


1 +d.z~ 1 1 + d,z -1 ’ * * 1 + d z" 1 

i c n 


(14) 


= (z) D 2 (z) . . . D n (z) 


(15) 


C l<z) c 2 (z)_ *„(.) 


E. (z) C (z) C . (z) 
l i n- t 


(16) 


The difference equations for Eqs. (14) and (16) are: 


C i (k) = a Q e.(k) - djC^k-1) 

C 2 (k) = Cj(k) + e 2 C 1 (k-l) -d 2 C 2 (k-l) 


e (k) = C (k) + e C (k-1) - d (k-1) 
n - l n n- 1 n 


( 17 ) 


18 

LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D225642 -II 



LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 


LMSC-HREC D225642-II 


A block diagram of this form is shown in Fig. 2-9. If each separate transfer 
function is programmed in the direct form and the output of one block is re- 
garded as the input to the next one. Then the number of storage registers 
for this form is 


Data 

Coefficients 

Multiplications 

Additions 

Transfers 


n 

m + n+ 1 
m + n+ 1 
m + n 
n 



Fig. 2-9 - General Block Diagram of the Cascade Programming Form 


2.5.4 Normal Programming 

Normal programming is based on the partial-fraction expansion of 
Eq. (8): 


D(z) = 


1 + dj z 


- 1 


i+d 2 z 


- 1 


+ . . . + 


_ n 

1 + d z 
n 


(18) 


20 


LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 




LMSC-HREC D225642-II 


where the are the poles of D(z) and the are the residues of D(z) at those 
poles. The difference equations for this form can be obtained by defining 


This leads to 


D(z) 


C . (z) C.(z) C (z) 

— = 4. — — x + II 

E. (z) E . (z) ' ' * E.(z) 

1 1 1 


Cj{k) = Pje.tk) - d^^k-1) 

c 2 ( k ) = P 2 e i < k ) " d 2 c 2^ k_1 ^ 


c (k) = p e.(k) - d c (k-1) 

n r n l n n 

n 

e o (k) = Z! c j (k) 

j = l 


(19) 


( 20 ) 


A block diagram of the normal form is shown in Fig. 2-10. The corresponding 
storage registers for this form are: 


Data 

n 

Coefficients 

2n 

Multiplications 

2n 

Additions 

2n - 1 

Transfers 

n 


th 

The equipment and operation requirements for implementing an m 
order polynomial in z in numerator and n order polynomial in z 1 in 
denominator of the function D(z) are listed in Table 2-2 for the various 
forms. The total number of storage registers is the sum of those re- 
quired for data, coefficients, multiplications, additions and data transfer 
storage. The total number of computer operations is the sum of multiplica- 
tion, addition and data transfer operations. Assuming each operation 
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Table 2-2 

EQUIPMENT AND OPERATION REQUIREMENTS 
FOR D(z) REALIZATION 



Forms 

Storage Registers 

Direct 

Canonical 

Cascade 

Normal 

Data 

m + n 

n 

n 

n 

Coefficients 

m + n+ 1 

m + n+ 1 

m + n + 1 

2n 

Multiplications 

m + n + 1 

m + n+ 1 

m + n + 1 

2n 

Additions 

m + n 

m + n 

m + n 

2n - 1 

Transfers 

m + n 

n 

n 

n 

Total Storage Registers 

5m + 5n+ 2 

3 m + 5n + 2 

3m + 5n+ 2 

8n - 1 

Total Computer 
Operations 

3m + 3n+ 1 

2m + 3n+ 1 

2m + 3n+ 1 

5n - 1 
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requires the same execution time, the computing time is proportional to the 
total number of computer operations. In general, each of these computer 
operations requires a different execution time. This table plus the fact that 
complex poles and zeros must be realized indicate that the canonical form of 
programming is desired for implementing the linear networks of the high 
frequency cutoff filter. 


2.6 QUANTIZATION ERROR 
2.6.1 Introduction 

Most numbers cannot be represented exactly in a digital system because 
of the finite word length of the system. The difference between the exact 
number and its representation in the digital system is referred to as the quan- 
tization error. This quantization error depends upon the type of mechanism 
used to represent the exact number. Generally, there are two types of quan- 
tization mechanisms used in digital devices — the truncation quantizer, and the 
roundoff quantizer. An input-output characteristic of a truncation quantizer 
is shown in Fig. 2-11. It represents only the first £ bits of an exact number, 
x. Therefore, the quantization error will not exceed 

h = + 2" £ (21) 

where h is called the quantization granularity. 

The input-output characteristics of a roundoff quantizer are shown in 
Fig. 2-12. For this quantizer, the maximum error due to quantization is 


h/2 = +2 1_£ 


( 22 ) 
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In order to compare quantization errors in digital systems programmed 
in canonical form and in direct form, expressions for the error are required. 
These are derived in the following section under the following assumptions: 

1. The digital filter is realized with fixed-point arithmetic 

2. Input signal to the filter is normalized such that the 
maximum input signal never exceeds unity 

3. All filter variables are scaled such that no overflows 
occur on addition. 

Since quantization is a nonlinear phenomenon, the present study on quan- 
tization is primarily concerned with obtaining an upper bound on the system 
output error. 


2.6.2 Quantization Error for Canonical Form 

Consider the digital control system as shown in Fig. 2-13 where D(z) is 
the pulse transfer function of a digital system (digital filter or digital con- 
troller). G, (s) is the transfer function of a zero-order hold and G(s) is the 
' ho 

transfer function of a plant. Assume that the digital filter has the following 
n 1 * 1 order pulse transfer function 


D(z) 


- 1 - 2 . . - m 

a Q + a l z + a 2 z + ... + a m z 

1 + b, z ^ + b,z ^ + . . . + b z n 

1 c n 


( 21 ) 


where m < n. If the digital system is programmed in the canonical form as 
shown in Fig. 2-8, and if each output of a time delay block is denoted as a 
state variable then the state equation can be represented as shown in Eq. (22). 
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x j (k+ 1 ) = x 2 (k) 
x 2 (k+l) = x 3 (k) 

x n _ 1 (k+i) = x n (k) 

n 

x n (k+l) = e.(k) - b n+1 _. Xj (k) (22) 

j = 1 

where x^(k) is the state variable x^ evaluated at time kT, and e(k) is the input 
evaluated at time kT. An additional equation specifying the output e Q (k) as a 
linear combination of the state variables is 


e (k) = a e. (k) + 
o' ' o 1' ' 


n 

£ 

j=i 


(a , , . - a 
n+l-j 


b . . )x.(k) 

o n+l-j' j x 


(23) 


Since 

e^k) = r(k) - c (k) (24) 

therefore Eqs. (22) and (23) can be rewritten as 

Xj (k+1) = x 2 (k) 
x 2 (k+l) = x 3 (k) 

x . (k+1) = x (k) 
n-1 n n 

■ x n (k+ 1 ) = r(k) - c(k) - b n+l-j x j^ k ) 

j=l 

n 

e o (k)=a o 0(k)-c(k)] + Yj (a n+l-j' a o b n+l-j )x j (k) (25) 

j = l 
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Furthermore, assume that the z-transform of the plant-hold combina- 
til 

tion has the following £ order transfer function 


a' + a’, z * + al, z * + . . . a r . z * 

TT , , _ o 1 2 1 

H( 2 ) ~ 

1 + b'j z + b' 2 z + . . . b^ z 

The state equations for Eq. (26) in the canonical form becomes 

y x (k+1 ) = y 2 (k) 
y 2 (k+l) = y 3 (k) 


(26) 


£ 

V r 1 -l< k ' 1) = e o <k) ■ 2 b j+l-j ^J <k) 

£ 

c(k) = a' o e Q (k) + £ (a- +1 _. - a^ b' + 1 _.) y.(k) 

j=l 


(27) 


Now the effect of input, output and multiplication quantization on system 
performance will be considered. Input quantization is caused by the conversion 
of input signal, e^(k), into a set of discrete levels. Output quantization is caused 
by the fact that not all the available bits of digital signal (e.g., e Q (k)) are con- 
verted to their analog equivalents. Multiplication quantization results when the 
product of two £ bit numbers is reduced to £ bits. 

The same control system with quantizers is shown in Fig. 2- 14, where 
the blocks labeled Q, Q^q» Qjj, are quantizers with quantization granu- 
larity hj, hji' • •• > respectively. The state equations for the digital 

system with quantizers become: 
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Xj(k+1) = x 2 (k) 
x 2 (k+l) = x^(k) 

x . (k+1) = x (k) 
n- 1 n 

x n (k+1)= ' Q ln [ b n x l (k) ] ' Q 1 (n- 1) [ b n- 1 X 2 (k) ] "• 

-°12 [ b 2 x „-l (k) ] - Q U [Vn (k >] +Q 10 [ r(k) ‘ '= lk) l 

e o <k)=Q 2n I a n x l (k> ] +Q 2(n-1) [ a n-l x 2 (k) ] + " + Q 2 1 [ a l x n (k) ) 

+ Q 20 | a o[- Q ln [ b n x l <k) ] Q 1 (n- 1) [ b „-I X 2 (k) l 

... - Q u |bj x n (k) j + Q 10 |r(k) - c(k)jj| (28) 


The bar over the state variables indicates that quantization effects have been 


included in the system. Q. b.x.(k) indicates the quantized value of b. x.(k). 

i 1 l 1 J * 1 J 

b^ x j(k)| may be defined as 


Q i [ b i X / k H = b i x j< k ) + ' ^ k ) 


where q^(k) is the error caused by the quantizer operating on b^x.(k) at the 
sampling instant t =kT. Similarly 


( 



Q. a 

Q j [ b 4 X m (k) ] 


x o 



= a oVm (k) + a o q j (k) + q i (k) 
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Expanding Eq. (28) yields 


x^k+1) = x 2 (k) 
x 2 (k+l) = x 3 (k) 


x . (k+1) = x (k) 
n- 1 n 


n 


c n (k+l) = r(k) - c <k) - ^ b n + 1 - j X j^ k ^ 

j=l 


n 


+ q io (k) ■ Z ^ 

j=i 


e o (k) = a Q [r(k) - c (k) ] + £ (a n+1-j ' a o b n+1 ..j> x.(k) 

j = l 


n 


n 


+ 2 q 2j (k> + q 20 (k) ‘ a o Z q lj (k) + a o q lO (k) 

j=l j=l 


(29) 


The quantized state equations for H(z) is 


yj(k+l) = y 2 (k) 


y 2 (k+l) = y 3 (k) 


X. 

^(k+D = Q [% (k) i ’ Z b i+l-j 7j( k ) 

l 

= %( k ) “ Z b kl-iM k) + q < k) 
j = l 
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c(k) = a^Q 


K< k >] + E - a oVi-j> ^ 
j=l 


= a' e 
o o 


(k) 


+ E (a i+l-j ' a o b «+l-j )y j (k) + a o q(k) 

j=l 


( 30 ) 


Since the effect of quantization on the system response is of interest, 
we may define error variables as 

y ^(k) = Xj, (k) - x^(k) 
u (k) = e Q (k) - e Q (k) 

(3 i ( k ) = y^k) - y i (k) 

v (k) = c (k) - c (k) (31) 

Subtracting Eq. (25) from Eq. (29) and Eq. (27) from Eq. (30) yields 


y x (k+ 1 ) = r 2 ( k) 
y 2 (k+ 1 ) = r 3 (k) 


y n _i(k+D = y n (k) 

n n 

y n (k+l) = - v(k) - 2 b n + 1 - j yj< k > + q 10 (k) ‘ H q lj 

j=l j =1 

n 

u(k) = -a o v(k) + (a n+1 _.- a o b n+1 _.)y.(k) 

j=l 

n n 

+ ]£ q 2j (k) + q 20 (k) ’ a o 2 q lj (k) + a o q 10 (k) 

j=l j=l 
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(3 ^k+1) = p 2 (k) 

P 2 (k+1) = (3 3 (k) 

P i-l (k+1 * = P^ k ) 

l 

P £ (k+1) = u(k)-]T b^ +1 _.3.(k) + q{k) 

j = l 
i. 

v(k) = a^u(k) + 0j ( k ) + a^q(k) (32) 

j=l 

Taking the z -transformation and assuming zero initial conditions, Eq. (32) may 
be expressed as 

zFjfz) = r 2 (z) 

zr 2 (z) = r 3 (z) 

zr „-l (z) = r n (z) 

»r n (z) =-vw - £ b n+ i_. riz) + q 10 ( 2 ) - ^Q,jU) 
j=l j=l 

n 

U<z) = -a o V(z)+ 2 (a n+l-j ' a o b n+l-j> 

j=l 

+ £ Q zj U) + Q 20 W - a o £ Q,.(z) + a o Q 10 (z) 

j=l j=l 
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z B^z) = B 2 (z) 
z B 2 (z) = B 3 (z) 

Z B *-l^ = 

l 

z B^z) = U(z) - £ b' B (z) + Q(z) 

j=l 

i 

V(z) = a Q U(z) + 2 (a' + 1 _. - a^b £+1 _.)B.(z) + a^Q(z) (33) 

j=l 


The steady state output error due to quantization can be obtained by 
simultaneously solving Eq. (33) for V(z) and then by applying the final value 
theorem (i.e. lim v(k) = lim (1-z *)V(z)) 
k — *■ oo z — *• 1 

For a = a' =0 
o o 


lim v(k) = 
k — * oo 




oo 



( l + S b i) lira f Z) q 2i lk) + q(k) ] | 

i=l k — *» oo l i=0 

( 1 + Z b .) ( 1 + E b i) + (£ 'i)(E h) 

i=l i=l i=l 


(34) 


Note that the maximum value of any of the q^(k) cannot exceed one quantization 
granularity k^. Therefore, the maximum steady state quantization error in 
the system output is: 
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lim v(k) 
k — fcco 


(X>i)|(i>i) (±*n) + ("±\){± h 2i + h )! 

= i =1 1=1 1=0 ' i=l \=0 

( 1+ t b i)( 1+ E b 0 + (E *.) (E*i) 

i=l i=l V i=l ' 'i=l ' 


(35) 


2.6.3 Quantization Error for Direct Form 

An upper bound on the steady state quantization error for the digital system 
programmed in the direct form can be obtained using a similar technique. A 
general block diagram of the direct programming form of a digital system whose 
transfer function is characterized by Eq. (21) was presented in Fig. 2-7. From 
Fig. 2-7, the state equations can be written as 


Xj(k+1) 


*0 1 0 ... 0 0 


x j (k) 


0 

X2(k+1) 


o 

o 

o 

o 


x 2 (k) 


0 

• 


• 

0 


• 


• 

• 

x . (k+1 ) 
n-1 


0 0 0 ... 0 1 


n . (k) 
n-1 


0 

x n (k+l) 

— 

000... 00 


x n (k) 

+ 

1 

*n+l <k+1 


010... 00 


*n+l (k > 


0 

x „+2 (k+1 > 


o 

o 

pH 

o 

o 


x n+2 (k) 


0 



0 ! 


• 


• 

x 2n-l (k+1) 


o 

o 

o 

o 

I—* 


x 2n-l< k > 


0 

x 2n ( k+1 ) 


000... 00 


X- (k) 

2n' ' 


0 


e.(k) 


0 

0 


0 

0 

0 

0 


u 1 a 


e 0 (k) 


( 36 ) 
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or 


x (k+1) = Ex(k) + F e.(k) + G e Q (k) 


(37) 


where E, F and G are defined by the equivalence of Eqs. (24) and (25). The 
output e Q (k) of the direct programming form is 


n 


e o (k) = 


i=l 


a. x . . (k) 
l (n+l)-i ' 


i=l 


V(2„ + l)-i< k)+a o a i (k) (38) 


Now the effect of signal and coefficient quantization on system performance 
will be considered. The digital system with quantizers is shown in Fig. 2-15. 

Let x(k) represent the state of the digital system with quantization included. 

Then the state equations of the system are 


x(k+l) = E x(k) + F Q I0 [e.(k)j + G e Q (k) 


(39) 


n 


n 


6 o (k) " °2i f a i X (n+l)-i (k) ^ ’ S Q li ^ b i X (2n+l)-i (k) ^ 

i=l i=l 


+ ^20 t a 0 ^10 [ e i^ k ^J 


Expanding Eq. (39) 

x(k+l) = E x(k) + F e.(k) + G e Q (k) + F q 1Q (k) 
n n 

a o (k) = E a i* (n+ i)-i 00 - E b i*<2n+l)-i (k)+a ° e i <k) 
n n 

+ £ q 2i (k) - 2 qii( k ) + a 0 q 10 (k) (40) 

i=0 i=l 
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Let 

y (k) = x(k) - x(k) . (41) 

and 

u(k) = e Q (k) - e Q (k) . 

From Eqs. (37), (38), (40) and (24), the following system error equations are 
obtained 


y(k+l) = 
u(k) = 

+ 


E y (k) + G u(k) + F q 1Q (k) - F v(k) 


(42) 


n 


n 


£ a i T(n+l)-i (k > - £ b i ’'(2n+l)-i< k) * a o v < k > 
n n 

Z q 2i ( k ) - Z q u (k) + a Q q 1Q (k) 

i=0 i=l 


The quantization error in the system output can be obtained by taking the z- 
transformation of Eq. (42) 


z £ (z) 
U(z) 


ET (z) + G U (z) + F Q in (z) - F V(z) (43) 

n 


— 10 ' 
n 


£ a i r (n +1 )-i( a )- S b. r (Zan)rt M -a o Y (a ) 


n n 

+ Z Q 2 i (x) " Z Q ii< z ) + a 0 Q 10 (z) 

i=0 i=l 
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Together with the corresponding error equations for H(z) (from Eq. (33)), the 
steady state output error due to quantization for the direct programming form 
can be obtained. For a Q = a^ = 0, the maximum steady state quantization error 
in the system output for the direct programming form is 


lim v(k) 
k — *> oo 



1 

— (44) 


Comparing Eq. (44) with Eq. (35) shows that the canonical programming 
form gives less quantization error than the direct form for those cases shown 
in Table 2-1 on page 10 if all the quantizers have the same quantization 
granularity. Note that h is a function of digital word length and the coeffi- 
cients a, b, a' and b' are functions of sampling rate and system time con- 
stants; therefore, Eqs. (35) and (44) may be used to estimate the digital word 
length and sampling rate required to achieve a given system specification. 

2.7 ANALYTICAL DESCRIPTION OF DIODE BRIDGE CIRCUIT* 


The nonlinear circuit consists of a diode bridge and an integrating ampli- 
fier with an input capacitance C equal to the feedback capacitance as shown 
in Fig. 2-16. 

2.7.1 Amplifier Transfer Function 

Let z be the input voltage to the capacitance C of the amplifier with 
grid voltage e and gain factor K, then 

o 


$ 

Acknowledgement is made to Dr. W. Trautwein of Lockheed -Huntsville for 
the analytical description of the diode bridge. 
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1 


(from lag filter) 



Fig. 2-16 - Nonlinear Part of Cutoff Filter 
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c ir( z - e s ) 


c w<*o- V 


and 

Substituting yields 

_d_ 

dt 

For K large, we obtain 


= -Ke. 


/ , e 0\ 

_d_ 


e 0 \ 

( 2+ k ) 

dt 

( e 0 + 

Tc 


e 0 « - 2 


i.e., the amplifier acts like an inverter. 

2.7.2 Diode Operation 

It is assumed that the bridge diodes have ideal switching characteristics 
that can be described by 

x , 


R = 


i oo if e < 0 


( 0 if e >0 

where e is the voltage across the diode. 

2.7.3 Diode Bridge Switching Logic 

The switching logic is governed by the respective values of the input 
voltages (ej,e 2) and the bridge output z. 

There exist six different conditions for the diodes which result in three 
different bridge outputs z as listed in Table 2-3. 
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Table 2-3 

BRIDGE OUTPUT AS A FUNCTION OF RESPECTIVE 
INPUT AND OUTPUT VOLTAGES 


Case 

No. 

Input -Output 
Voltage Relation 

Resulting Diode States 

Resulting 
Bridge Voltages 

Resulting 
Bridge Output 

Diode No. 

1 

2 

3 

4 

X 

y 

la 

e 2 > e l > z 

N 

C 

C 

N 

e l 

e 2 

z = e j 

lb 

z > e i > e 2 

C 

N 

N 

C 

e 2 

e l 

z = e j 

2a 

e i > e 2 

c 

N 

C 

N 

e 2 

e l 

z=e 2 

2b 

z > > e i 

N 

C 

N 

C 

e l 

e 2 

z = e 2 

3a 

e i > z > e 2 

C 

N 

N 

N 

e 2 

e l 

z = z(0) 

3b 

e 2 >z >e { 

N 

C 

N 

N 

e l 

e 2 

o' 

N 

II 

N 


C = conducting; N = nonconducting 


Note that both Diodes 3 and 4 are blocked in Case 3 which makes it impossible 
for the capacitor to discharge. This results in a "hold" -type bridge output in 
Case 3. The value z(0) in Table 2-3 is to be taken at the most recent switching 
time of the diodes. 


The complete bridge plus capacitor plus amplifier signal transfer can 
now be described by the relationship: 



if 

e 2 

> e l > z(0) 

or 

z(0) 

- e l > e 2 

if 

e l 

> e 2 > z(0) 

or 

z(0) 

>e 2 > e 

if 

e l 

> z(0) > e 2 

or 

e 2 

> z(0) > 


( 45 ) 
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Equation (45) can be readily programmed for a digital computer by a number 
of logical decisions ; e. g. , 

A MX = A MAXI (el, e2) 

AMN = A MINI (el, e2) 

If (Z ‘ GT • A MX) eO = - A MX 
If (Z • LT • AMN) eO = - AMN 
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Section 3 

HIGH FREQUENCY CUTOFF FILTER PERFORMANCE 

3.1 ANALOG FILTER PERFORMANCE 

The analog realization of the cutoff filter is analyzed to determine its 
input-output characteristics. This analysis included the dc input plus high 
frequency signal. 

3.1.1 DC Signal Only 

An analog simulation of this filter is shown in Fig. 3-1. The primary 
reason for conducting this experiment was to test the dc characteristics of a 
circuit with an amplifier using capacitive input and feedback elements. As one 
may see from the test results, the filter performs at dc quite satisfactorily. 
The various results are shown in Figs. 3-2 and 3-3. 

TEST PROCEDURE 

The analog circuit of Fig. 3-1 was simulated using the parameter values 
shown in Table 3-1. The dc input-output characteristics were found by two 
methods, slope input and step input. 

Table 3-1 

PARAMETER VALUES USED IN THE SIMULATION OF FIG. 3-1 
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;o: 


L© — ^ — I 


+ 100 V 


+ 100 V 




Pot Values 


Notes: 


Value 

No. 

Value 

Wj 2 / 10000 

7 

K/10 

2^/100 

8 

.0035 

Wl 2 /l0000 

9 

.0035 

U 2 /100 

10 

v L /<ioo + v L > 

W 2 /100 

<2 

11 

Vj/000 + VJ 


1. Time scale : l/l 

2. K = 

3. Vl is absolute value of voltage at which 
A-Z is to be limited to prevent overload 
during testing. 

4. Numbers preceding integrators and 
amplifiers indicate gain settings for 
the integrators and amplifiers. 


Fig. 3-1 - Analog Simulation Diagram for HFCF Filter 
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Time 

(sec) 


Fig. 3-3 - Step Function Response 
of the Cutoff Filter for 
Filter Parameters Shown 
in Table 3-1 


Fig. 3-2 - Input-Output Characteristics 
of the Cutoff Filter 
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Slope Input Test 

The input to the circuit was made to follow a ramp function from -25 volts 
to +25 volts in approximately 520 seconds. The output, after an initial transient, 
was also a ramp function, slightly delayed, and of almost equal magnitude. It 
was thus ascertained that the dc operation of the filter was approximately linear . 


Step Input Test 

To construct an input-output curve of the filter (shown in Fig. 3-2), a 
series of step function inputs was made to the filter. These inputs were varied 
from -25 volts to +25 volts in 2.5 volt intervals. The outputs were obtained 
and the steady-state output values were recorded by a digital voltmeter print- 
out. These input-output voltages are the points shown in Fig. 3-2. As may be 
seen, a constant dc offset of -0.05 volts was encountered. The step function 
response is shown in Fig. 3-3 for an input of 20 volts. The transient response 
shows that the filter clips oscillations in the same way as in previous high- 
frequency tests (Ref. 1). The particular rise time, settling time, etc., of the 
step response in Fig. 3-3 is of course a function of the parameters chosen and 
given in Table 3-1. The step response could obviously be altered by varying 
these parameters within the bounds given in Fig. 2-1. 

3.1.2 DC Plus High Frequency Signal 

In the previous section, the results of dc input-output characteristics 
of a high-frequency cut-off filter were presented. These characteristics indi- 
cated a linear relationship between input and output dc levels. The question 
arises, will this relationship remain linear if the filter is subjected to an input 
consisting of dc plus a high frequency sinusoid whose amplitude is greater 
than the dc level. The answer was determined to be negative. In order to 
maintain the desired linear relationship, it is necessary that the dc level be 
greater in amplitude than the negative peak value of the time varying input. 

This is easily accomplished, however, as indicated by the block diagram of 
Fig. 3-4, which illustrates a modification of the cut-off filter which allows the 
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dc transfer function to be held essentially constant (as shown in Fig. 3-5). Some 
error was introduced to the data plotted in Fig. 3-5 by the necessity to read values 
from strip chart recordings, which are reproduced in Fig. 3-6. The requirement 
that the dc bias signal be chosen to maintain a positive waveform into the original 
filter is not restrictive and very much simplifies a describing function approach 
to analysis since the number of possible input cases is then reduced by one 
half. 

3.2 DIGITAL FILTER PERFORMANCE 

A flow chart of the digital program for the high frequency cutoff filter is 
shown in Fig. 3-7. The linear filter networks are programmed in the canonical 
form. 


The program has the capability to read in sampling period, cutoff frequency, 
filter frequencies and damping ratios. It then calculates the corresponding HFCF 
coefficients, frequency characteristics of both continuous filter and digitized 
filter, and time responses for different inputs. Both the frequency response and 
time response subroutines of the program contain plot routines. 

3.2.1 Transient Response 

Figures 3-8 through 3-17 show the transient response of the digital high 
frequency cutoff filter for various combinations of ^ j and t, The purpose 
of this is to study the effect of damping ratio upon the transient response of the 
filter. The figures show that the amount of overshoot depends on the damping 
ratios C, ^ and t, If one of the t, is critically damped, i.e., damping ratio 

equal to one, there is no or little overshoot. If both £, ^ and t, ^ are underdamped, 
i.e., damping ratios less than one, overshoot occurs. The peak overshoot is 
the first overshoot. The output also oscillates around the final value. The over- 
damped system is slow-acting but does not oscillate about the final value. The 
underdamped system reaches the final value faster than the overdamped case, 
but the response oscillates about the final value. The amount of permissible 
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Fig. 3-7 - Flow Chart of Digital Program to Calculate HFCF Coefficients, 
Frequency Characteristics and Time Responses 
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Fig. 3-8 - Time Response of Digital Filter 
(Step Input = 20 at t > 0) 
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Fig. 3-9 - Time Response of Digital Filter (Step input = 20 
at t >0) 
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Fig. 3-10 - Time Response of Digital Filter (Step input = 20 
at t ? 0) 
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Fig. 3-11 - Time Response of Digital Filter (Step input = 20 
at t >0) 
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Fig. 3-12 - Time Response of Digital Filter (Step input = 20 
at t ? 0) 
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Fig. 3-13 - Time Response of Digital Filter (Step input = 20 
at t >0) 
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Fig. 3-14 - Time Response of Digital Filter 
(Step Input = 20 at t > 0) 
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Fig. 3-15 - Time Response of Digital Filter 
(Step Input = 20 at t > 0) 
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Fig. 3-16 - Time Response of Digital Filter 
(Step Input = 20 at t > 0) 
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Fig. 3-17 - Time Response of Digital Filter (Step Input = 20 
at t > 0) 
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overshoot determines the desirable values of damping ratios. Figure 3-18 shows 
Ihe responses of the HFCF for a constant C, 2 of 0.7 and three different values 
°f !» j* -A family of log magnitude curves and phase curves for these combina- 
tions of damping ratios are also plotted in Figs. 3-19 and 3-20. These curves 
show the effects of damping ratio upon the frequency response of the linear 
portions of the digital filter for input frequencies equal to or less than one half 
the sampling frequency. 

3.2.2 DC Plus High Frequency Signals 


A number of simulations were run for various inputs to investigate the time 
response and input-output relations of the digital high frequency cutoff filter. 

Some time responses of the filter are shown in Figs. 3-21 through 3-24. The 
following filter parameters were used. 


HFCF Parameters 


Cutoff Frequency, Cl> c 

3 rad/sec 

Sampling Frequency, f 

s 

25 sample s/i 

“l 

4.5 rad/sec 

tl 

0.5 

W 2 

6.0 rad/sec 

^ 2 

0.6 


Figures 3 22 and 3-23 show that for dc plus a single high frequency signal, 
the dc signal is transmitted and the high frequency signal is blocked by the digital 
filter. 


A comparison of filter response to a step input (Fig. 3-21) with the response 
of the filter to a step input plus high frequency signal (Figs. 3-22 and 3-23) shows 
that the high frequency signal has caused a decrease in rise time and an increase 
in overshoot. The step response of the filter in the presence of one, two, three 
and five tones of high frequency signals is shown in Fig. 3-24. Figure 3-25 
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Fig. 3-18 - Time Response of High Frequency Cutoff Filter 
for Three Different Values of 
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Fig. 3-20 - Phase Diagram for the HFCF for Several Combinations of £ and / 

1 2 
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Fig. 3-21 - Time Response of High Frequency Cutoff Filter 
(Step Input = 20 at t<^0) 
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Fig. 3 22 - Time Response of High Frequency Cutoff Filter (Input = 20 + 20 sin 12t) 
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Fig. 3-23- Time Response of High Frequency Cutoff Filter (Input = 20 + 20 sin 18t) 
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shows the dc input-output characteristic of the digital high frequency cutoff 
filter. Note the linear relationship between dc input and output level. 

3.2.3 DC Plus Low Frequency Signal 

Figure 3-26 shows the response of the digital filter to a dc plus a low 
frequency signal for the following filter parameters 


Cutoff Frequency, ( 0 c 
Sampling Frequency, f 

s 



3 rad/sec 
25 samples/sec 
4.5 rad /sec 
1.0 

6.0 rad/sec 
0.7 


As may be observed, both the dc and the low frequency signals are trans- 
mitted by the digital filter. Hysteresis effect occurred because of the high 
frequency cutoff characteristics of the filter logic. 


3.2.4 Low Frequency Signal 

Figures 3-27, 3-28 and 3-29 show the time responses of the digital filter 
to sinusoidal signals of frequencies 0.4 C 0 c> 0.6 <J 0 c and 0.8 co c , respectively. 

The filter output is characterized by hysteresis. This effect is greater 
as the signal frequency approaches the cutoff frequency. The hysteresis is due 
to the high frequency cutoff characteristic of the digital logic. Note that the 
steady state low frequency magnitude characteristic of the HFCF can be described 
by the hysteresis loop as shown in the sketch on the following page. 
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max 

HFCF Hysteresis Characteristic 

3.2.5 Low Frequency and Multiple High Frequency Inputs 

Figure 3-30 shows the time history of the filter output in response to 
an input containing multiple high frequency components in addition to a low 
frequency component. Figure 3-31 is a plot of the filter output and the low 
frequency component of the input. The filter parameters are the same as 
those used in Section 3.2.3. We see from these plots that the low frequency 
component is passed in a "quasi-sampled" manner. 

It has been shown (Ref. 1) that due to the high frequency attenuation 
pattern, the filter output has an appearance similar to that resulting from 
the superposition of a single high frequency signal on that of a low frequency. 
This single high frequency is the lowest expected frequency above u > c . 

Hence, regardless of the high frequency content of the input, analytical studies 
to describe the nonlinear low frequency effects due to quasi-sampling can be 
based on a filter input consisting of a single high frequency signal superimposed 
on a signal of low frequency. 
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Fig. 3-26 - Time Response of High Frequency Cutoff Filter 
(Input = 20 + 20 sin 1.2t) 
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Fig. 3-27 - Time Response of High Frequency Cutoff Filter 
(Input = 20 sin 1.2t) 
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Fig. 3-28 - Time Response of High Frequency Cutoff Filter 
(Input =20 sin 1.8t) 
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Fig. 3-29 - Time Response of High Frequency Cutoff Filter (Input =20 sin 2.4t) 

77 

LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 




Filter Input and Output 


LMSC-HREC D225642 -II 



Time (sec) 

Fig. 3-30 - Time Response of High Frequency Cutoff Filter 
(Input = 20 (sin 0.6t + sin6t + sin 12t + sin I8t) 


78 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D225642 -II 



Time (sec) 


Fig. 3-31 - Time Response of High Frequency Cutoff Filter 
(Input = 20 (sin 0.6t + sin6t + sin 12t + sinl8t) 
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3.2.6 High Frequency Signal 

Figures 3-32, 3-33 and 3-34 show the time responses of the digital 
filter to sinusoidal signals of frequencies ut^, 4o» c and 8<J c> respectively. 

The following filter parameters were used: 

U) c 3 rad/sec 

f 25 samples/sec 

lo j 4. 5 rad/ sec 

0.5 

a >2 6.0 rad/sec 

( 2 °' 6 

These figures show that when a high frequency signal alone appears at the 
filter input, sharp cutoff occurs and the input is subdued. 

Figures 3-35 and 3-36 show that the time response of the digital filter 
output to two out-of-phase high frequency signals is periodic. The filter 
parameters used were 

u 2 7 r rad/sec 

c 

< * 37 r rad/sec 

K x o - 5 

47T rad/sec 

*2 °’ 6 

For a sampling period of 0.04 second, the period of the output wave- 

form is 5 seconds and for a sampling period of 0.08 second, the period of the 
output waveform is 10 seconds. These are equal to the least common multiple 
of the periods of the high frequency signals and the sampling period. Note 
also that the magnitude of the output is very small compared to the magnitude 
of the input. 
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Fig. 3-32 - Filter Response to High Frequency Signal 
(Input = sin 3t, = 1) 
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Fig. 3-33 - Filter Response to High Frequency Signal 
(Input = sin 12t, =4) 
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Fig. 3-34 - Filter Response to High Frequency Signal 
(Input = sin24t, t * ) c = 8) 
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Fig. 3-35 - Filter Response to Out-of-Phase High Frequency Signals 
(Input = sin 367Tt + l4cos87Tt, Sampling Period = 0.04 sec) 
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Fig. 3-36 - Filter Response to Out-of-Phase High Frequency Signals 

(Input = sin 3.6 7rt + l4cos87rt, Sampling Period = 0.08 sec) 
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Section 4 

CONCLUSIONS AND RECOMMENDATIONS 

The problem of digital implementation of a cutoff filter developed by 
Lockheed Missiles & Space Company's Huntsville Research & Engineering 
Center is approached with consideration to word length, sampling rate, ac- 
curacy requirements, computing time and hardware restrictions. 

Upper bounds for the steady state system output error due to quantization 
for digital control systems containing a digital network programmed both in the 
direct form and in the canonical form were derived. This was accomplished by 
defining a set of error equations in the z domain then applying the final 
value theorem to the solution. Quantization error was found to depend upon the 
digital word length, sampling rate, and system time constants. The error 
bound developed may be used to estimate the digital word length and sampling 
rate required to achieve a given system specification. 

Computing time and hardware requirements for four possible programming 
forms for the linear portions of the filter were determined. From the quantiza- 
tion error accumulation, computing time and hardware point of view, and the fact 
that complex poles and zeros must be realized, the canonical form of program- 
ming seems preferable. 

Finally, the basic input-output characteristics of the digitized high frequency 
cutoff filter were investigated. The filter demonstrated its ability to pass signals 
with minimum distortion in a given frequency band while minimizing sig na ls whose 
frequencies lie outside the band. 

Since the most important performance criterion of all is overall system 
stability, it is recommended that this study be continued so that the following 
stability problems can be investigated: 
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1. Development of a general method (e. g. , Lyapunov's 2nd method) 
to analyze the stability of a given control system containing the 
high frequency cutoff filter. 

2. Extension of the above developed method to analyze the stability 
of sampled-data control systems containing the digitized high 
frequency cutoff filter. 

3. Investigation of sampling frequency effects on the stability of a 
sampled-data system containing the digitized filter. 


It is also recommended that the analytical results on quantization error 
be verified by simulating the high frequency cutoff filter using fixed point 
arithmetic. 
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Appendix A 

HYSTERESIS OF DIGITAL HIGH FREQUENCY 
CUTOFF FILTER 
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Appendix A 

In the case where no high frequency signals are present the passage of 
a low frequency signal is characterized by hysteresis. To describe the phase 
and gain characteristics due to hysteresis, the describing function of the non- 
linearity is derived. From Fig. 3-29 it is observed that the digital output e Q 
is neither an even nor odd function and therefore must be described by both 
the sin n co t and cos n 6) t terms of the Fourier series. Assuming e has a 
period of Zir and the sampling rate is n, then the Fourier series of the out- 
put is given by 


e o (k) = 


a /2 
o' 


00 

Ek 

i=i l 1 


cos 


/ Zir 
[n+l 


k )‘ 


+ b. cos 

l 


( 2 TT 

\ n+l 


k H 


where k = 0, 1, . . . , n . 


The coefficients of the first H harmonics are given by Ref. 8. 


n 


a. 

l 


n+l 


E 

k=0 


e (k) cos 
o' ' 


2?rk 

n+l 


b. 

l 


9 

n+l 


eJk) sin i 

k=0 ° n+1 


For input = E 1 sin t, the describing function of the nonlinearity is then 
given by 



A - 1 -CL/ 
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whe re 


Zj = (a* + b* ) * 

and 

-1 a l 
<f) , = tan rr— 

b l 


Using the above method, the describing function for several values of co 
is determined and listed in Table A-l for the following digital filter parameters. 


Cutoff Frequency, a> c 

2jt rad/sec 

Sampling Frequency, f 

s 

25 samples/sec 

“l 

1.5 co 

c 

tl 

0.5 

0) 2 

2.0 co 

c 

^ 2 

0.6 


Table A-l 

VALUES FOR DESCRIBING FUNCTION 


a )/< o c 

E 1 

a l 

b l 

Z 1 

*1 

.2 

10 

-.121 

9.86 

9.86 

42 ' 

.3 

10 

-.251 

9.67 

9.67 

1 ° 30 ' 

.4 

10 

-.425 

9.16 

9.18 

2 ° 38 ' 

.5 

10 

-.564 

8.61 

8.63 

3 ° 45 ' 

.6 

10 

00 
1 — i 

• 

1 

7.80 

7.85 

5 ° 16 ' 

.7 

10 

-.708 

6.36 

6.40 

6 ° 21 ' 

.8 

10 

-.738 

4.75 

4.80 

00 

o 

oo 

.9 

10 

-.512 

2.54 

2.60 

11 ° 32 ' 


A-2 
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